Composite boson mapping for lattice boson systems.

We present a canonical mapping transforming physical boson operators into quadratic products of cluster composite bosons that preserves matrix elements of operators when a physical constraint is enforced. We map the 2D lattice Bose-Hubbard Hamiltonian into 2×2 composite bosons and solve it within a generalized Hartree-Bogoliubov approximation. The resulting Mott insulator-superfluid phase diagram reproduces well quantum Monte Carlo results. The Higgs boson behavior in the superfluid phase along the unit density line is unraveled and in remarkable agreement with experiments. Results for the properties of the ground and excited states are competitive with other state-of-the-art approaches, but at a fraction of their computational cost. The composite boson mapping here introduced can be readily applied to frustrated many-body systems where most methodologies face significant hurdles.